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I. INTRODUCTION 



The subject of this work is the extended three-electron Hylleraas integrals involving ^/rfj 
and l/rf terms, which appear in matrix elements of relativistic operators in the Breit-Pauli 
Hamiltonian. These integrals are defined as 

J An J Att J An 

^ni-l 712-1 ns-l 714-1 rtg-l na-l /i \ 

'23 '31 '12 '1 '2 '3 5 \^) 

with one of Ui equal to —1, and all other are nonnegative integers. They have been 

l|,0,3 andby Yan et al. Q. 



studied in detail in a series of papers by King 



Their approach 



is based on the expansion of l/r^" in an infinite series of some orthogonal polynomials. 
The resulting multiple summation is performed with the help of convergence accelerators. 
In some special cases |ll this expansion can be avoided and Eq. (jT)) can be expressed in 
terms of the two-electron Hylleraas integral. Using these methods leading relativistic and 
QED corrections to energies ^,0], isotope-shift 7|, and ^f-factors ^ have been calculated 
to a high degree of accuracy. The analytic approach to Hylleraas integral has been so far 
.u* ,e. .u..e.M. P.O.. and HO, ,„ |^ we. aMe .o oMa. a„ ana,^,c exp.e.,o„ to. 
a more general integral with exponents in r^, as well as in rij. However, their expression 
is quite complicated, involves multivalued dilogarithmic functions, and thus is of limited 
use. In parallel Remiddi in obtained a simple analytic expression for the Hylleraas 
integral with rij = 0. Recently, together with Remiddi we derived recursion relations jll| for 
Hylleraas integrals for arbitrary large powers of and rjj. This result allows for a convenient 
calculation of the nonrelativistic wave function of the lithium atom and light lithium-like 
ions. 

In this work we present an analytic approach which allows for a fast and high precision 
calculation of extended 3-electron Hylleraas integrals. This sets the ground for improvin g th e 
precision of theoretical energy levels by including higher order relativistic corrections ■ 
We obtain here the closed form expression for the master integral /(— 1, 0, 0, 0, 0, 0) and 
derive recursion relations for increasing values of arguments rii of /(— 1, n2, ns, n4, ns, ne). 
The other integrals, such us /(ni, ?t,2, 'n.3, — 1, ris, ng), can be obtained by using already de- 
rived recursion relations for /(ni, n2, Tis, n4, ns, ng) with nonnegative followed by a one 
dimensional numerical integration with respect to the corresponding parameter Wi. 
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II. RECURSION RELATIONS FOR rr^ INTEGRAL 



Our derivation is based on integration by parts identities which are commonly used for 



the calculations of multi-loop Feynman diagrams We follow here the former work 
and first consider the integral G 



G(mi,m2,m3;m4,m5,m6) = J d^ki J (fk2 j d^k^{k\ + u\) "^^{kl + ul 

{kl + ul)--^ {kl, + wl)-^^ {kl, + wl)--^^ {kl, + wl)--%2) 

which is related to / by: /(O, 0, 0, 0, 0, 0, 0) = (7(1, 1, 1, 1, 1, l)|ui=„2=n3=o- The following 9 
integration by part identities are valid because the integral of the derivative of a function 
vanishing at infinity vanishes, 

= id(i, j) = / £k^ [ d^k2 [ d^ks^lkj {kf + ul)-^ 

J J J dki 

{kl + uir' {kl + ui)-\ki, + wir' {kl, + wir' {ki, + wir'] , (3) 



where z,j = 1,2,3. The reduction of the scalar products from the numerator leads to the 
identities for the linear combination of the G-function. If any of the arguments is equal to 
0, then G becomes a known two-electron Hylleraas type integral, Eq. ()B6|) . The explicit 
form of all 9 identities is presented in Eq. (jAlll and the whole derivation presented below is 
performed with the help of MATHEMATICA program for symbolic computation. 

In the first step of deriving recursion relations we take the difference id(3, 2) — id(2, 2) 
and use it as an equation for G(l, 2, 1, 1, 1, 1), 

G(l, 2, 1, 1, 1, 1) {ul -ul + wl) = G(l, 1,1,0, 1, 2) - G(l, 1, 1, 1, 0, 2) 

+G(1,1,1,1,1,1)-G(1,2,0,1,1,1) 
+G(1, 2, 1, 0, 1, 1) - 2 G(l, 1, 1, 2, 1, 1) wl 
+G{l,lXl,l,2){wl-wl-wl). (4) 

Similarly, the difference id(2, 3) — id(3, 3) is used to obtain G(l, 1, 2, 1, 1, 1), 

G(l, 1, 2, 1, 1, 1) {ul -ul- wl) = G(l, 0, 2, 1, 1, 1) - G(l, 1, 1, 0, 2, 1) 

-G(l, 1,1,1, 1,1) + G(1, 1,1, 1,2,0) 
-G(1,1,2,0,1,1) + 2G(1,1,1,2,1,1)«;? 

+G{l,l,l,l,2,l) {wl + wl-wl). (5) 



These two equations are used now to derive recursions in 712 and 713. With the help of the 
trivial formula 

roo p-ur 1 

3, (6) 



/ du 

Jo r 



one integrates with respect to ui, which lowers the first argument rii to —1. Next, one 
differentiates with respect to U2,U3,Wi,W2, W3 at U2 — — to generate arbitrary powers 
of ri3, ri2, Ti, r2, ra and obtains quite long recursion relations for n2 and na, 

1 



/(-I, n2 + 2, na, 714, 715, 7i6) 



wa 



/(-1, 7i2, na + 2, 7i4, ris, 7i6) 



7l4 (^4 - 1) (7l2 + 1) /(-I, ^2, na, 7l4 " 2, 7l5, 7l6 + 1) 
-7l5 (715 - 1) (7l2 + 1) /(-I, ^2, 7l3, 7l4, - 2, 7l6 + 1) 

+n6 (1 +n2 + 2rz4 + uq) (n2 + 1) /(-I, n2, ns, 724, ns, ne - 1) 
-na (na - 1) /(-I, n2 + 2, 7ia - 2, 714, 715, 7i6 - 1) 
+7i4 (714 - 1) Tie /(-I, n2 + 2, 7ia, n4 - 2, 715, - 1) 

-2 7Z4 (712 + 1) Wi /(-I, 77.2, 77.3, n,i - 1,775,776 + 1) 

-2 7i6 (n2 + 1) wi /(-I, n2, na, n4 + 1, 715, 7i6 - 1) 
-2n4neWif{-l,n2 + 2,7ia,n4 - I,n5,7i6 - 1) 

+776 ^i^l /(-I, n2 + 2, 773, n4, U^, Hq - 1) 

+2 n5 (772 + 1) W2 /(-I, n2, na, 774, 775 - 1, 775 + 1) 

-(2 + 772 + 2 774 + 2776) (n2 + 1) 7^3 /(-I, 772, na, 774, ns, 775) 

+ (773 - 1)773703 /(-I, 772 + 2,n3 - 2,774,775,776) 

-774 (774 - 1) W3 /(-I, 2 + 772, n3, -2 + 774, 775, 775) 
+2 (772 + 1) wi W3 /(-I, n2, na, 774 + 1, n^, uq) wi W3 

+2 774 Wi W3 /(-I, n2 + 2, 773, 774 - 1, 775, 776) 

+ (772 + 1) (Wi -wl + wl) /(-I, n2, n3, 774, n5, 776 + 1) 

+5^4 W3 r(77a + 775 - 1, 772 + n6 + 1, -1; «^2, W3, 0) 

-5„3 wa r(-l + 774 + 775, 776, n2; Wl + W2, W3, 0) 

-(772 + n6 + 1) 5„4 r(773 + 775 - 1, 772 + TIq, -1] W2, W3, 0) 
+776 Sna ^{^A + n5 - 1, 776 - 1, n2; Wi + W2, Ws, 0) 

+(772 + 1) r(n6 - 1, na + 774 - 1, 772; W3, Wi, 0) , (7) 

-772 (n2 - 1) ns /(-I, -n2 - 2, 773 + 2, 774, 775 - 1, 776) 



W1W2 



+ (?23 + 1) ^4 (n4 - 1) /(-I, n2, na, n4 - 2, ns + 1, ne) 
+ (n3 + 1) ris (1 + + 2 724 + ng) /(-I, ^2, ng, 7^4, rig - 1, rig) 
-(ng + 1) ne (rig - 1) /(-I, ris, ^3, ^4, + 1, rig - 2) 
+^4 (n4 - 1) ^5 /(-I, ^2, ris + 2, n4 - 2, ^5 - 1, tiq) 
-2 (na + 1) Hi wi /(-I, n2, ns, 714 - 1, + 1, ng) 
-2 (ns + l)n5 wi /(-I, n2, ns, n4 + 1, ns - 1, uq) 
-2 n4 ns wi /(-I, n2, 2 + ns, n4 - 1, ns - 1, uq) 
+ (ns + 1) {wl + wl- wj) /(-I, n2, ns, n4, ns + 1, tiq) 
+^5 wl /(-I, n2, ns + 2, ^4, ns - 1, ng) 
+n2 (n2 - 1) W2 /(-I, n2 - 2, ns + 2, n4, ns, Uq) 
-(ns + 1) (ns + 2 n4 + 2 ns + 2) w;2 /(-I, ^^2, ^4, '^s, ^e) 
-n4 (n4 - 1) W2 /(-I, n2, 2 + ns, n4 - 2, ns, ng) 
+2 (ns + 1) wi W2 /(-I, n2, ns, n4 + 1, ns, uq) 
+2n4 wi W2f{-l,n2,n3 + 2,n4 - l,ns,n6) 
+2 (ns + 1) ne W3 /(-I, n2, n4, ns + 1, ng - 1) 
W2 r(ns + ns + 1, n2 + ng - 1, -1; W2, W3, 0) 
-5„2 W2 r(n4 + ng - 1, ns, ns; Wi + W3, W2, 0) 
+ (ns + l)Sne + n4 - 1, ns - 1, ns; wi, W2, 0) 
-(1 + ns + ns) r(ns + ns, n2 + ng - 1, -1; W2, W3, 0) 
+ns 5n2 r(n4 + ng - 1, ns, ns - 1; Wi + W3, W2, 0) , (8) 

where (5„ denotes Kronecker delta 5n,o and F is a two-electron Hylleraas integral, which 
is defined in Eq. (jBljl . These recursions assume that the values of /(— 1, 0, 0, n4, ns, ng), 
/(— 1, 1, 0, n4, ns, ng), /(— 1, 0, 1, n4, ns, ng) and /(— 1, 1, 1, n4, ns, ng) are known. We cal- 
culate master integrals for the last three cases explicitly and express them in terms of 
two-electron Hylleraas integrals as in 

;(_l.,l,0,0.0)^/^/^/''^'' 



wl {wl - wfj W3 wf 



4 71 r|s ri r2 rs 
\n^ = ^T{0,0,-l;w2,W3,0), (9) 
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Ti f d^r2 f e 



1 



■wi ri—W2 r2—W3 ra 



4Wi W2 

-2Li2 
1 



An 



W2 



W2 + W3 
W2 



rh ri r2 
21 ^2 



+ 2Li2 1 



1^1 + ^2 + 1^3 
W3 

W2 



+ 2U2 



- 2 Lio 1 - 



W2 



W2 + W3 
Wi + W3 



Wi + W2 + W3/ \ W2J V W2 

r(-l, 0, -1; ws, W2, 0) - r(-l, 0, -1; wi + W3, W2, 0) 



(10) 



^3 r d^ r2 f d^ e 



„_,i,o,o.o.o,^/^^/fL^/ 



■wi T-1-UI2 r2— U13 r3 



4 vr r|3 ri2 n r2 



r(-l, 0, -1; W2, W3, 0) - r(-l, 0, -l;wi + W2, W3, 0) 



(11) 



where Li2 is a dilogarithmic function, Eq. ()C13|) . The recursion relations in 77,4,77.5,716 are 
obtained by differentiation with respect to wi, W2 and W3 



/(-I, 772 , 713, 714,715, Tie) = — 



-7I4 (7I4 - 1) /(-I, 7l2, 7I3, 7I4 - 2, 7I5, Uq) 
+2 7I4 Wi /(-I, 7l2, 713, 714 - 1, 715, Tie) 

r(Tl5 + TI3 - 1, Tie + TI2 - 1, -1; W2, W3, 0) 
-5„3 r(Tl4 + TI5 - 1, Tie, -1; Wi + W2, W3, 0) 
-Sn2 r(Tl4 + Tie - 1, 775, -1; Wi + W3, W2, O) , (12) 

for (T12,T13) G {(1,1),(1,0),(0,1)}. 

The calculation of the integral /(—I, 0, 0, ri4, TI5, Tie) is much more elaborate and we have 
to return to integration by parts identities, see Appendix A. These are 9 equations, which 
we solve against the following Xi=i^g unknowns at U2 = U3 = 0, 





= G(l,2,l,l,l, 


l)Ui, 


X2 


= ^(1, 1,2,1,1, 


l)ui, 


X3 


= G(l,l,l,l,2, 


1)mi, 


X4 


= ^(1, 1,1, 1,1, 


2)ni, 


X5 


= G(l,2, 1,1,1, 


1), 


Xe 


= ^(1, 1,2, 1,1, 


1), 


Xr 


= ^(1, 1,1, 1,2, 


1), 


Xs 


= G(l,l,l,l,l, 


2), 


X9 


= G(2, 1,1, 1,1,1). 
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(13) 



The solution for X7 and Xs is 



^(1, 1,1, 1,2,1) = -^(2, 1,1, 1,1,1) 



2wi 



^(1, 1,1, 2, 1,1) 



wf + W2 



Wn 



2wl 



4wf W2 



(14) 



^(1, 1,1, 1,1, 2) 



2 9 

-G(2, 1, 1, 1, 1, 1) 7^ - G(l, 1, 1, 2, 1, 1) 



2wi 



+^(1, 1,1, 1,1,1) 



3wj — W2 + 



4 wf 



4wf W3' 



(15) 



where 



G'(l,l,l,0,2,l)u;2_^(^^^^^^^^2,0)m 
-G(2, 0, 1, 1, 1, 1) wl - G(2, 1, 1, 1, 1, 0) wl 
-G(l, 1, 1, 0, 1, 2) wl + G(l, 1, 1, 1, 0, 2) wl 
-G{2, 1, 0, 1, 1, 1) wl + G(2, 1, 1, 1, 0, 1) wl 



(16) 

We now use exphcit form of two-electron integrals in Eq. ()B6jl and integrate both equations 
with respect to mi, 

\f +{2wl + wl - wl) /(-I, 0, 0, 0, 0, 0) 

2WiW2^ 

-Wl {wl + wl - wl) /(-I, 0, 0, 1, 0, 0)1 , (17) 



/(-I, 0,0, 0,1,0) 



/(-I, 0,0,0,0,1) 



-F + (2 wl - wl + wl) /(-I, 0, 0, 0, 0, 0) 



2wj W3 

-w, {wl - wl + wl) /(-I, 0, 0, 1, 0, 0) 



(18) 



where 
F 



dui F{ui) 



2 Lio - 



W2 



+U2 



Wl 
W2 



\W2 + W3 



Li2^1 
- Li2 



W2 

Ws 



Lis 1 



+Li,(l-^ 

W2 



Lis 1 



W2 + W3 

Wl + W3 



+ Lis 



Wl + W2 

W3 
W2 



WiJ 



W1+W2+ W3 



Lis 



W3 



W1+W2 + W3 



+ lnf^)ln 

\W3 



W1+W2+ W3 



W2 / \W3/ \ W2 + W3 

r(0, -1, -2; W2, Wl + W3, 0) - r(0, -1, -2; ^3, Wi + W2, 0) 
+r(0, -1, -2, 0; Wl, W2) - r(0, -1, -2, 0; Wi, W3) 
+W2 r(-l, 0, -1; Wl, 0, W2) - W3 r(-l, 0, -1; Wl, 0, W3) 
+W2 r(-l, 0, -1; 0, W2, W3) - W3 r(-l, 0, -1; 0, W3, W2) . 



(19) 



Next, we multiply both equations by powers of Wi to eliminate any Wi from the denominator 
and differentiate with respect to wi, W2 and w^. This leads to the following recursion relations 
in and rig of the /-function 



1 

2Wi W2 

—2 n4 (n4 — 1) W2 /(—I, 0, 0, 77,4 — 2, 77,5 + 1, ng) 



/(-I, 0, 0, ^4, ^5 + 1, ne) = K {n4 - 1) (^4 + 2 ^5) /(-I, 0, 0, - 2, ^5, ng) 



-714 (3774 + 4775 + 1) /(-I, 0, 0,774 " 1, ^5, ^^g) 
+4774^1 77;2 /(-I, 0,0,774 - 1, 775 + Ij^-e) 
+ (774 + 1) (775 - 1)775 /(-I, 0,0, 774, 775 - 2,776) 
-2 (774 + 1) 775 W2 /(-I, 0, 0, 774, 775 - 1, 776) 
-(774 + 1) (776 - 1) ri6 /(-I, 0, 0, 774, 775, 776 " 2) 
+2 (774 + 1) 776 W3 /(-I, 0, 0, 774, 7I5, 776 - 1) 

+[(3 774 + 2 775 + 2) + (774 + 1) 77;2 - (774 + 1) W3] 

X/(-l,0,0,774,775,776) 

-(775 - 1) 77.5 Wi /(-I, 0, 0, 774 + 1, 775 - 2, 776) 

+2 775 wi W2 /(-I, 0, 0, 774 + 1, 775 - 1, 775) 

+ (776 - 1) rie Wi /(-I, 0, 0, 774 + 1, 775, 776 - 2) 
-2 776 Wi W3 /(-I, 0, 0, 774 + 1, 775, 776 - 1) 
-Wi {wf + wl- wl) /(-I, 0, 0, 774 + 1, 775, 776) 



+F(774,775,776)J, (20) 
1 

2wi W3 

2 (774 - 1) 774 77;3 /(-I, 0, 0, 774 - 2, 775, 776 + 1) 



/(-I, 0,0, 774, 775, 776 + 1) = 2 ... r74 (774 - 1) (774 + 2 776) /(-I, 0,0, 774 - 2,775,776) 



-774 (3 774 + 4 776 + 1) Wl /(-I, 0, 0, 774 - 1, 775, 776) 

+4774 77;iW3/(-l,0,0,774 - 1,775,776 + 1) 

-(774 + 1) (775 - 1) 775 /(-I, 0, 0, 774, 775 - 2, 776) 

+2 (774 + 1) 775 W2 /(-I, 0, 0, 774, 775 - 1, 776) 

+ (774 + 1) (776 - 1) n-Q /(-I, 0, 0, 774, 775, 776 - 2) 

-2 (774 + 1) 776 W3 /(-I, 0, 0, 774, 775, 776 - 1) 

+ [(3 774 + 2 776 + 2) - (^4 + 1) 7^2 + (^4 + 1) ^^3] 



X/(-l, 0,0,714,715,^6) 

+ (77,5 - 1) 7l5 Wi /(-I, 0, 0, 7Z4 + 1, 7Z5 - 2, 715) 
-2 7l5 Wi W2 /(-I, 0, 0, 7l4 + 1, 7l5 - 1, 7l6) 

-(7i6 - 1) ne wi /(-I, 0, 0, 7i4 + 1, n5, uq - 2) 

+27l6WiU'3/(-l,0, 0,7l4 + I,7l5,7l6 " 1) 
-Wi {wl -wl + wl) /(-I, 0, 0, 7l4 + 1, ^5, 716) 
-F(7l4,n5,7l6) , (21) 



where 



F{n,,n^,n^) = {-d^,Y^ {-8^,^^ {-d^.Y^ F 

= (7l6 - ns) r(-l, 7l5 - 1, 7l6 - 1; 0, W2, W3) 
+W2 r(-l, 7l5, 7l6 " 1; 0, W2, W3) 

-W3 6n^ r(-l, riQ, 7Z5 - 1; 0, W3, W2) 

-(715 - 1) Sne r(7l4 - 1, 0, 715 - 2; Wi, 0, ^2) 
+W2 Sne r(7l4 - 1, 0, 715 - l;Wi, 0, ^2) 

+ (7i6 - 1) 5„s r(7i4 - 1, 0, Tie - 2; Wl, 0, W3) 

-W3 r(-l + 7l4, 0, 7l6 - 1; Wl, 0, ^3) 
+r(7l5, 7l4 + 7l6 - 1, -2; W2, + W3, 0) 

-r(7l6, 7l4 + 7l5 - 1, -2; W3, Wl + W2, 0). (22) 

What remains is the calculation of /(—I, 0, 0, n4, 0, 0). In the following we derive a differen- 
tial equation for h{wi) = /(—I, 0, 0, 0, 0, 0), from which we obtain /(—I, 0, 0, 714, 0, 0). The 
solutions for 

(7(2, 1,1, 1,1, 1)^X9, (23) 

and for the difference 

Xi«r'-^5 = 0, (24) 

form two algebraic equations, which however are too long to be written here. They in- 
volve the terms ^(2, 1, 1, 1, 1, 1), G{2, 1, 1, 1, 1, l)ul G{1, 1, 1, 1, 1, 1), ^(1, 1, 1, 1, 1, 1) u^^, 
G(l, 1, 1, 2, 1, 1), G(l, 1, 1, 2, 1, 1) M^^^, and the known two-electron terms, where one of the 
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arguments of G-function is equal to 0. We integrate both equations in ui from e to oo, 
approach the hmit e — > and drop In e 



/oo 
dm G{2, 1,1, 1,1,1) 

/oo 
duiulG{2, 1,1, 1,1,1) 

/oo 
du,G{l,l,l,l,l,l) = h{w,), 

POO 

du, ^(1, 1, 1, 1, 1, 1) = giw,) - /(I, 0, 0, 0, 0, 0), 



2 ' 
hjwi) 
2 ' 



h'{w^) 



fOO 

du,G{l,l,l,2,l,l) = - 

r 1 .-^r(^ 1 1 9 1 n 3'iw,) 1 0,0, 0,0,0) 
I % G{1, 1,1,2, 1, 1) = ^ + — , (25) 



where 



h{w,) = /(-I, 0,0, 0,0,0), 

= [^e-^^^^-^^"^^^{lnr23 + j)r^^'r^2'r^'r^'r^\ (26) 



and 



/(I, 0,0, 0,0,0) = ^In 



Wi {Wi +W2 + W3) 



(27) 



.{wi + W2) {wi + W3). 

The set of both equations forms first order differential equations for h{wi) and g{wi). We 
eliminate g{wi) to obtain the following second order differential equation for h{wi) 

w wl h"{wi) + ivi [4 wi {wi — wi — wi) + w] h'iwi) 

+ + 4wi {wl -wl- wl) - w] h{wi) = R{wi), (28) 

where 

'Wl + w^ 



R(wi) = Wl W2 ln( 1 + — ) + Wl W3 ln( 1 + — ) + (wl - wl) ln( - 

V ^2/ V W3J \Wi + W2 

21 f Wi{Wi+W2 + W3) \ 2 2^ TP 

+2wi In ^ — + (wo - ^3) F, 



(29) 



and 



4, 4, 4 o22 o22 o22 
W = Wl + W2 + W^ — 2WiW2 — 2W2W3 — ZWiW^ 

= -(-Wl + ^2 + W3)(wi -W2 + W3)(wi +W2- W3)(wi + W2 + W3). (30) 
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Two linearly independent solutions of the homogenous equation are: 



hi{wi) 



1 



{Wi +W2- W3) {Wi - W2 + W3) 

{-Wl + W2 + W3) {Wi + W2+ W3) 
4^2 W3 



Wl y/W2 W3 

where is a complete elliptic integral of the first kind as defined in Eq. 
Wronskian W is 

W = hi{wi) h'2{wi) - h[{wi) h2{wi) 



(31) 



11), and the 



2n 



W Wl 



(32) 



where w is defined in Eq. ()3()|1 . The solution in Eq. is valid for wi in the range 

\w2 — W3\ < Wl < W2 + W3, because the elliptic integral K has a branch cut for arguments 
exceeding 1. We use the identity in Eq. ()C8|) . to obtain solution hi of the homogenous 
equation for wi > W2 + W3 



hi{wi) 

Wl y {wi + W2- Ws) (Wi -W2 + W3) 

and h2 for wi < \w2 — w^l 

o 

h2{Wi) 



K 



4 W2 W3 



{Wi +W2- W3){Wi -W2 + W3) 



(33) 



K 



Wl J {-Wl + ^2 + W3) {Wi +W2+ W3) 



4:W2 W3 



{-Wl +W2 + W3){wi +W2+ W3) 



(34) 



The solution of the inhomogeneous equation is obtained by Euler's method of variation of 
the constant, 



h{wi) 



hi{wi) f^2+w3 ^R{w')h2{w') h2{wi) 



27r 



dw 



+ 



W' 



27r 



dw 



,R{w') hi{w') 



\w2-W-i\ 



W' 



(35) 



There is no additional term being a solution of homogenous equation, because h{wi) is 
finite for all values of Wi, but not hi{wi) and h2{wi). Therefore this is the right solution. 
Having obtained /(— 1, 0, 0, 0, 0, 0) = h{wi) and /(— 1, 0, 0, 1, 0, 0) = —h'{wi) we calculate 
/(—I, 0, 0, n4, 0, 0) = h{wi,n4^) = {—dwi)^^h{wi) recursively. The inhomogeneous differential 
equation is differentiated n-times with respect to Wi, to obtain 

h{wi,n + 2) = -^-^\-{n-3){n-2)^{n-l)nh{wi,n-4:) 
WWi ^ 

+ {n — 2) {n — 1) n (13 — 17n + 6 n^) Wi h{wi, n — 3) 

-{n -l)n [(14 - 25 n + 15 n^) wl-2{n- if w J h{wi, n - 2) 

+2nwi {3 — hn + \{\v?')w\ + {—1 + 3n — An'^)ws h{wi,n — l) 
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+ 



(4 + 10 n + 15 n^) + (1 - n^) + 2 (1 + 3 n + 6 n^) Ws h{wi,n) 



+Wi 



(5 + 6 n) + (1 + 2 n) - 2 (3 + 4 n) h{wi,n + 1) 



+R{wi, n 



where = w| + W3, = w| — w1 and 



(36) 



R{wi, n) 



R{wi,2) 



4(n-3)! 5w2{n-2)\ 



4(r2-3)! 5^3 (n- 2)! 



n-2 



{Wi + W2Y ^ {Wi + W2Y ^ {Wi + ^ 

4(n-3)! 2{w2 + w:if{n-l)\ A{w2 + w^) {n - 2)\ 



{Wi + Wis)" 2 

4(n-3)! 
iwl-wl)F{n,0,0), 



(wi +W2 + 



{Wi +W2 + W^Y ^ 



^ (41^1-^3^) {n-l)\ ^ {Awl - ^i) {n - 1)! 



for n > 2, 



R{wi), 

-{W2 + W3) + 



Awl 



Wn 



- 2 



-W3 ln( 1 + — 

5W2 ^ 



Wi + W2 

— Awi In 

5^3 



+ 



Awl 



Wo 



{W2 + WsY) 
Wi + W2 + W3 Wi + W3 
Wi {Wi + W2+ W3) 



W2 In ( 1 + — 

V W2 



{Wi + W2) {Wi + W3) 
2 {W2 + W3Y A {W2 + W3) 



+ {wl-wl)F{l,0,0), 



Wi + W2 Wi + W3 {Wi +W2 + W3Y W1+W2+ W3 



Awl 



w 



3 ^ 4w;i 



Wn 



{W1 + W2Y {wi + w3y 

+ {wl-wl)F{2,0,0). 



+ 4 In 



Wi {Wi +W2 + W3) 
{wi + W2) (wi + W3) 



(37) 



In the case Wi 



w. 



smg 



W2 + W3 or \w2 — Wsl, the recursion in Eq. (jHBj) is not numerically 



stable. Therefore, one instead of this recursion, calculates the recursion exactly at Wi 
what corresponds to setting w = in Eq. 

h(wi, n + 1) = - 



wi \{5 + 6n)wf + {1 + 2 n)w^- 2 {3 + A n) wj Ws 
{-(n - 3) (n - 2)^ (n - 1) n h{wi, n - A) 

+ {n — 2) {n — 1) n (13 — 17n + 6 n^) Wi h{wi, n — 3) 

-{n-l)n {lA-25n + 15n'^)wf-2{n-lYws h{wi,n-2) 

+2nwi (3 — 5n + lOn^) wl + {—I + 3n — An^) Ws h(wi,n — l) 



+ 



-(4 + 10n + 15n^)w^ + {l-n^)wl + 2{l + 3n + 6 n^) w^ w^ h{wi,n) 



+R{wi, n 



(38) 
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where wi = Wsmg- This completes the recursion relations for the extended Hylleraas integral 
with r^3^. 

III. NUMERICAL EVALUATION 

We pass now to numerical implementation of recursions and integration of the master 
integral in Eq. (j35p . All the computation is performed with extended precision arithmetic, 
namely quad and sextuple precision. Even higher precision, the octuple one, is used for 
checking numerical accuracy. The starting point is the master integral. One needs to calcu- 
late it with the highest possible accuracy, because the recursions depend most significantly 
on the value of initial terms. The integrand in Eq. ()35p is a product of the function R and 
the complete elliptic integral K. The function R defined in Eq. has singularities only 
at Wi = and Wi = oo, and the complete elliptic integral has logarithmic singularities at 
Wi = Wsing = \u!2 — Ws\ or W2 + Ws which correspond to zeros of w in Eq. (jHOI)- In the 
following we assume that 1^2 — ^3 7^ 0. When \ w2 — W3\+ei < wi < W2 + W3 — €2, the integral 
in Eq. ()35|) can be performed by the Gauss- Legendre quadrature (l^ . We have verified that 
for €1 ~ 0.2 the integration with 100 points gives the accuracy of at least 32 digits if not 
more. For the cases Wi > W2 + W3 + €2 and Wi < \w2 — W3I — ei, the integration contour 
is deformed on the complex plane to avoid the singularity. So this contour consists of two 
lines on a real axis, joined by a half-circle with origin at the singular point and integration 
is performed independently on each part using the Gauss-Legendre quadrature. When wi 
is close to the singular point, we first obtain h{wsing ± e), next we calculate h{wsing) by 
matching the Taylor expansion from recursion in Eq. ()38|) with h{wsmg ± e) and in the last 
step we again use recursion in Eq. to obtain h{wi). For W2 — it's = we separate R in 
Eq. (j29|) into the part which is free of \\iWi and the part which is proportional to Inwi. The 
first part is integrated using the Gauss-Legendre quadrature an the second one is integrated 
using the Gauss quadrature adapted for the logarithmic weight function. Several numerical 
results for some selected Wi are presented in Table I. 

Considering recursions, all but one involve denominators limited from below. Only that 
in Eq. (j36|) for increasing ^4 has a denominator which can be arbitrarily close to zero. 
Therefore, if Wi ~ w^ing-, one instead of recursion in Eq. (j36p . uses the recursion in Eq. 
(PH|l and calculates h{wi,n) from Taylor expansion at wi = Wsmg- All other recursions are 
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TABLE I: Values of the master integral for selected wi,W2, and 



Wi 


W2 


W3 


/(-I, 0,0, 0,0,0) 


4.0 


1.0 


0.5 


1.243 735 828 073 620 173 310 981 564 244[-l] 


4.0 


1.0 


1.0 


9.855 133 136 060 504 470 218 647 797 889[-2] 


4.0 


1.0 


1.5 


8.181 412 007 841 597 436 460 514 476 518[-2] 


4.0 


1.0 


2.0 


6.983 588 391 604 680 181 982 031 823 035 [-2] 


4.0 


1.0 


2.5 


6.077 218 287 692 100 226 048 417 176 715[-2] 


4.0 


1.0 


3.0 


5.365 400 720 042 544 709 716 176 264 176[-2] 


4.0 


1.0 


3.5 


4.791 010 346 652 078 406 517 300 908 585[-2] 


4.0 


1.0 


4.0 


4.317 729 831 064 450 749 511 048 756 748[-2] 


4.0 


1.0 


4.5 


3.921 185 585 221 614 693 378 573 393 156[-2] 


4.0 


1.0 


5.0 


3.584 332 630 993 527 980 351 431 968 712 [-2] 


4.0 


1.0 


5.5 


3.294 856 745 699 432 037 984 459 599 008[-2] 



calculated directly as in Eqs. (0 |Hl [121 EOl 1211) • They involve two-electron integrals F. The 
calculation of T including singular cases has recently been described in detail in Refs. 3| 
and jr3], and it does not pose any problem. Finally, several numerical results for three- 
eleven integ.^ .nvolvng powe. of and r, a.e presented . Table 11. Fo. comparison 
with the former results obtained in Ref. yj, we have taken the same m and Wi as in Table 

n 

II of this Ref. Our results agree up to the precision achieved in Ref. p|| with the one 
correction. In the fifth position instead of /(2, 1, 1, 3, 3, —2, 4.338, 4.338, 7.384), it should be 
/(1, 1, 2, 3, 3, -2, 4.338, 4.338, 7.384). 

Considering the extended Hylleraas integral with rf ^, we calculate it by numerical inte- 
gration with respect to Wi 

POO 

/(ni,n2,n3, -I,n5,n6) = / c?wi /(ni, na, ng, 0, ng, ne). (39) 



The recursion relations for /(^i, n2, ns, n4, ns, ng) with nonnegative Ui have been derived 
previously and they seem to be stable enough to perform this integration numerically. 
Since, we have not found in the literature the method of integration, which is adapted to two 
weight functions: the constant and the logarithmic one, we use the standard Gauss-Legendre 
quadrature. In Table III we present several numerical results. It is observed the significant 
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TABLE II: Three-electron Hylleraas integral involving l/r^2- ^^'^ 'ire from Table II of Ref. 
P]. Function / from this Ref. should be divided by (47r)'^ for comparison with our function /. 



Wl W2 W3 



f{ 


-1 


2 


2, 


2, 


2, 


2) 


2.700 


2.700 


2.700 


3.622 072 193 238 069 065 841 911 460 566 


-3] 


f{ 


-1 


2 


4, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


2.044 941 897 990 188 175 637 070 889 313 


-1] 


f{ 


-1 





2, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


8.560 152 684 198 427 372 519 849 562 718 


-3] 


f{ 


-1 





0, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


7.695 548 443 927 856 456 193 296 733 495 


-3] 


f{ 


-1 


4 


4, 


3, 


2, 


2) 


7.384 


4.338 


4.338 


2.516 457 130 304 929 175 434 829 560 592 


-6] 


f( 


-1 





0, 


2, 


2, 


2) 


3.000 


2.000 


1.000 


7.759 319 533 814 226 728 190 558 692 235 


-3] 


f( 


-1 





0, 


1, 


3, 


2) 


3.000 


1.000 


2.000 


1.528 428 874 506 937 507 531 543 743 291 


-2] 


/( 


-1 





4, 


2, 


4, 


3) 


2.000 


3.000 


4.000 


6.208 037 315 282 433 323 108 Oil 184 899 


-3] 


/( 


-1 


2 


2, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


4.036 629 272 285 446 411 970 138 933 470 


-2] 


/( 


-1 


2 


2, 


1, 


1, 


1) 


2.500 


2.500 


0.600 


1.025 702 855 657 754 018 359 340 659 240 


-1] 


/( 


-1 


2 


2, 


0, 


0, 


0) 


2.700 


2.900 


0.650 


3.674 068 373 009 625 515 617 159 197 784 


-2] 


f( 


1 

— 1 


9 


9 


1 

1 , 


9 




1 nnn 


1 nnn 


1 nnn 

1 .uuu 


A ^7fi 9Q^ Af\'i zL^I QSzL ^1 4 0"^ ^ nQ7 S7Q zL1 1 
^.0 t u zyo ^uo ^01 yo^ 01^ aoo uc/ 1 / y ^11 


91 


/( 


-1 


2 


2, 


2, 


2, 


2) 


1.000 


1.000 


1.000 


5.436 536 048 634 697 021 325 813 246 683 


2] 


/( 


-1 


2 


4, 


4, 


3, 


1) 


3.000 


2.000 


1.000 


6.126 463 692 215 932 446 059 888 955 061 


0] 


/( 


-1 


2 


4, 


1, 


0, 


0) 


1.000 


1.000 


1.000 


4.219 398 540 932 754 336 898 663 066 822 


2] 


/( 


-1 


4 


4, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


3.826 213 635 276 544 192 395 399 453 200 


0] 


/( 


-1 


4 


6, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


4.206 326 264 336 604 338 380 540 655 410 


1] 


/( 


-1 


6 


6, 


1, 


1, 


1) 


2.700 


2.900 


0.650 


1.886 948 258 407 236 970 462 772 961 418 


3] 



loss of precision due to the presence of Inwi/wi at large wi asymptotic. Therefore, precise 
integration requires in some cases, the subtraction of these terms. 

IV. SUMMARY 

An analytic approach is presented for the calculation of three-electron Hylleraas inte- 
grals involving one inverse quadratic power of inter-particle coordinate. This approach is 
based on exact recursion relations in powers of coordinates. These recursions involve initial 
terms and two-electron integrals. For the initial term /(—I, 0, 0, 0, 0, 0) as a function of Wi, 
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TABLE III: Three-electron Hylleraas integral involving 1/rf at wi = 2, tt)2 = 3, t^a = 4 



fff] n n — 1 






f(ci 1 n —1 




1 '^7fiQ8^9fi^ ^07 0*^9 1fi4f-9l 


f('n 9 0—1 




(\ Q49 9fiQ *^fiQ DQ"^ 71 9 1 0i^f 

U.i7T:^ ^Uc/ OUc/ UjU I XZ( ±UtJ|-OJ 






AO*^ 99*^ ^"^1 R^ ^ RQ^ 1 1 R[ 

(J.'4:UO ZZO ■4:01 OltJ OaO llO -OJ 


ffT) 4 0—1 






/(O, 5, 0,-1 


0,0) 


8.587 459 945 883 427 557[-3] 


/(O, 6, 0,-1 


0,0) 


1.598 496 287 482 975 980[-2] 


/(O, 7, 0,-1 


0,0) 


3.698 745 219 132 481 190[-2] 


/(O, 8, 0,-1 


0,0) 


1.036 442 843 454 920 448[-l] 


/(O, 9, 0,-1 


0,0) 


3.434 721508 609 856189[-1] 



one constructs a linear second order differential equation. Its solution is expressed as a 
one-dimensional integral over dilogarithmic and elliptic function which can be obtained 
numerically with arbitrary high precision. The two-electron Hylleraas integrals have already 
been derived in the literature and they also can be obtained with arbitrary high precision. 
These extended Hylleraas integrals are necessary for the calculation of relativistic effects 



in lithium and light lithium-like ions p, |6|, lD, |8||. One interesting goal is the high precision 
calculation of the lithium hyperfine splitting [18|, which can serve as a benchmark result 
for other less accurate methods. Moreover, it has recently become possible to derive nu- 
clear parameters such as the charge radius from the measurement of the isotope shift 
The hyperfine splitting is sensitive to the distribution of magnetic moment. Therefore the 
measurement of hfs in various odd isotopes of Li or light lithium-like ions may lead to the 
determination of the so called magnetic radius, which is very difficult to access experimen- 
tally. Even more interesting is the possible extension of this analytic method to beryllium 
and beryllium-like ions, the 4-electron systems. The use of Hylleraas basis set will allow for 
a high precision calculation of the wave function and, for example, various transition rates 
which are of astrophysical relevance. 



16 



V. ACKNOWLEDGMENTS 



We are grateful to Vladimir Korobov for his source code of the fast multi-precision arith- 
metics and the dilogarithmic function. We wish to thank Krzysztof Meissner for help in 
solving differential equation and Frederick King for presenting us his numerical results for 
some selected values of extended Hylleraas integral and for useful comments. This work was 
supported by EU grant HPRI-CT-2001-50034. 



[1] F.W. King, Phys. Rev. A 44, 7108 (1991); F.W. King, K.J. Dykema and A.D. Lund, Phys. 
Rev. A 46, 5406 (1992). 

[2] I. Porras and F.W. King, Phys. Rev. A 49, 1637 (1994); P.J. Pelzl and F.W. King, Phys. Rev. 

E 57, 7268 (1998); F.W. King, Adv. At. Mol. Opt. Phys. 40, 57 (1999). 
[3] P.J. Pelzl, G.J. Smethells, and F.W. King, Phys. Rev. E 65, 036707 (2002). 
[4] Z.-C. Yan, G.W.F. Drake, J. Phys. B 30, 4723 (1997); Z.-C. Yan, J. Phys. B 33, 2437 (2000). 
[5] F.W. King, J. Mol. Struct. (Theochem) 400, 7 (1997); F.W. King, et al, Phys. Rev. A 58, 

3597 (1998). 

[6] G.W.F. Drake and Z.-C. Yan, Phys. Rev. A 46, 2378 (1992); Z.-C. Yan and G. W. F. Drake, 
Phys. Rev. A 52, 3711 (1995); Phys. Rev. Lett. 81, 774 (1998); Phys. Rev. A 61, 022504 
(2000); Phys. Rev. A 66, 042504 (2002); Phys. Rev. Lett. 91, 113004 (2003). 

[7] Z.-C. Yan and G. W. F. Drake, Phys. Rev. A 66, 042504 (2002). 

[8] Z.-C. Yan, J. Phys. B 35, 1885 (2002). 

[9] D. M. Fromm and R. N. Hill, Phys. Rev. A 36, 1013 (1987). 
[10] E. Remiddi, Phys. Rev. A 44, 5492 (1991). 

[11] K. Pachucki, M. Puchalski and E. Remiddi, Phys. Rev. A. 70, 032502 (2004). 

[12] K. Pachucki, Phys. Rev. A to be published. 

[13] F.V. Tkachov, Phys. Lett. B 100, 65 (1981); K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. 
B 192, 159 (1981). 

[14] W. H. Press, S.A. Teukolsky, W.T. Vetterling and B.F. Flannery, Numerical Recipes in FOR- 
TRAN 77, 2nd ed. (Cambridge University Press, Cambridge, 1992). 
[15] R.A. Sack, C.C.J Roothaan and W. Kolos, J. Math. Phys. 8, 1093 (1967). 



17 



[16] V.I. Korobov, J. Phys. B 35, 1959 (2002). 

[17] F.E Harris, A.M. Frolov and V.S. Smith, Jr., J. Chem. Phys 121, 6323 (2004). 
[18] K. Pachucki, Phys. Rev. A 66, 062501 (2002). 

[19] M. Abramowitz and I. A. Stegun (1974), Handbook of Mathematical Functions (Dover, New 
York, 1964). 

APPENDIX A: INTEGRATION BY PARTS IDENTITIES 

The complete set of recursion relations as obtained from integration by parts identities 
is presented below. Function G is defined in Eq. (0) and id(z, j) = for i,j = 1, 2, 3, 



id(l,l) 


= -^(0, 


1,1, 


1,1,2) 


- ^(0, 1, 1, 1, 2, 1) + ^(1, 0, 1, 1, 1, 2) + ^(1, 1, 0, 1, 2, 1) 




-G{1, 


1,1, 


1,1,1) 


+ 2 G(2, 1, 1, 1, 1, 1) ul + G(l, 1, 1, 1, 2, 1) {ul -ul + wl) 




+G(1, 


1,1, 


1,1,2) 




id(2,l) 


= -^(0, 


1,1, 


1,1,2) 


- G(0, 1, 1, 1, 2, 1) + G(l, 0, 1, 1, 1, 2) + G(l, 1, 0, 1, 2, 1) 






1,1, 


0,2,1) 


+ ^(1, 1, 1, 1, 2, 0) - G{2, 0, 1, 1, 1, 1) + G{2, 1, 1, 1, 1, 0) 




+G(1, 


1,1, 


1,1,2) 


{ul -ul- wl) + G(2, 1, 1, 1, 1, 1) {ul + ul- wl) 




+G(1, 


1,1, 


1,2,1) 


1 2 2,2 2\ 
Ml — Mo + — o , 


id(3, 1) 


= -^(0, 


1,1, 


1,1,2) 


- ^(0, 1, 1, 1, 2, 1) + ^(1, 0, 1, 1, 1, 2) + ^(1, 1, 0, 1, 2, 1) 




-G(l, 


1,1, 


0,1,2) 


+ G(l, 1, 1, 1, 0, 2) - G(2, 1, 0, 1, 1, 1) + G(2, 1, 1, 1, 0, 1) 




+G(1, 


1,1, 


1,2,1) 


{ul -ul- wl) + G{2, 1, 1, 1, 1, 1) {ul ^u\- wl) 




+G(1, 


1,1, 


1,1,2) 


1 2 2,2 2\ 


id(2,2) 


= G(0,1 


1,1 


,1,2)- 


- G(l, 0, 1, 1, 1, 2) - G(l, 0, 1, 2, 1, 1) + G(l, 1, 0, 2, 1, 1) 




-G(l, 


1,1, 


1,1,1) 


+ 2 G(l, 2, 1, 1, 1, 1) ul + G(l, 1, 1, 2, 1, 1) {ul -ul^ w\) 




+^(1, 


1,1, 


1,1,2) 


{-ul^ul^w\), 


id(l,2) 


= ^(0,1 


1,1 


,1,2)- 


- ^(0, 2, 1, 1, 1, 1) - ^(1, 0, 1, 1, 1, 2) - ^(1, 0, 1, 2, 1, 1) 




+G(1, 


1,0, 


2,1,1) 


- G(l, 1, 1, 2, 0, 1) + G(l, 1, 1, 2, 1, 0) + G(l, 2, 1, 1, 1, 0) 




+G(1, 


1,1, 


1,1,2) 


{-ul ^ul- wl) + G(l, 2, 1, 1, 1, 1) {ul + ul- wl) 






1,1, 


2,1,1) 


{ul-ul + wl-wD, 


id(3,2) 


= ^(0, 1,1, 1,1,2) - 


- G{1, 0, 1, 1, 1, 2) - ^(1, 0, 1, 2, 1, 1) + ^(1, 1, 0, 2, 1, 1) 
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+G(1, 1, 1, 0, 1, 2) - Gil, 1, 1, 1, 0, 2) - Gil, 2, 0, 1, 1, 1) + 2, 1, 0, 1, 1) 

+G(1, 1, 1, 2, 1, 1) [ul -ul- wl) + G{1, 2, 1, 1, 1, 1) {ul + ul- wl) 

+G{1, 1, 1, 1, 1, 2) {-u\ + ul-wl + wl), 
id(3, 3) = G(0, 1, 1, 1, 2, 1) + G(l, 0, 1, 2, 1, 1) - ^(1, 1, 0, 1, 2, 1) - ^(1, 1, 0, 2, 1, 1) 

-G{1, 1, 1, 1, 1, 1) + 2 G{1, 1, 2, 1, 1, 1) ul + G{1, 1, 1, 2, 1, 1) i-ul + ul + wl) 

+G{l,l,l,l,2,l){-u\ + ul + wl), 
id(2, 3) = G(0, 1, 1, 1, 2, 1) + G{1, 0, 1, 2, 1, 1) - G{1, 0, 2, 1, 1, 1) - G{1, 1, 0, 1, 2, 1) 

-G(l, 1, 0, 2, 1, 1) + G(l, 1, 1, 0, 2, 1) - G{1, 1,1,1, 2, 0) + ^(1, 1, 2, 0, 1, 1) 

+G(1, 1, 1, 2, 1, 1) (-^2 + ul- wl) + ^(1, 1, 2, 1, 1, 1) (u^ + ul- wl) 

+G{1, 1, 1, 1, 2, 1) {-u\ + ul-wl + wl), 
id(l, 3) = ^(0, 1, 1, 1, 2, 1) - G(0, 1, 2, 1, 1, 1) + G{1, 0, 1, 2, 1, 1) - G{1, 1, 0, 1, 2, 1) 
1, 0, 2, 1, 1) + G(l, 1, 1, 2, 0, 1) - G{1, 1, 1, 2, 1, 0) + G{1, 1, 2, 1, 0, 1) 

+G(1, 1, 1, 1, 2, 1) {-u\ + ul- wl) + G(l, 1, 2, 1, 1, 1) [ul + ul- wl) 

+G(1, 1, 1, 2, 1, 1) i-ul + ul- wl + wl). 

(Al) 



APPENDIX B: TWO-ELECTRON INTEGRALS 

The two-electron integral F is defined by 

r(ni,n2,n3;ai,«2,«3) =1^1 ^ e""^^^""^^^"'^^^- r^"' r^-Vrr^ (Bl) 

J 4 TT J 4 TT 

In the simplest case of rii = ?t,2 = '^3 = it is 

r(0, 0, 0, ; «i, «2, «3) = , ^ w w ^ (B2) 

(ai + 0:2) (a2 + as) (0:3 + ax) 

In the general case of > 

, ,/d\"i/d\"2/d\"3 1 

r(ni,n2,n3;ai,a2,a3) 



dcti/ V dQ;2/ V dQ;3/ (oi + Q;2)(a2 + a3)(«3 + «i) ' 

(B3) 



and recursions relations for its evaluations have been derived in Refs. |15i|. The two-electron 
integral F for < can be obtained by the integration over Oj. Typical examples are 

^ / d\"2/ d ^"3 in(a, + a,)-ln(ai + a3) 
r(-l,ra2,n3;ai,a2,a3) = -75 — [--, — 7 ^ ^ , (B4) 

V da2/ V dQ:3/ Q!2 - Q!3 "2 + "3 
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r(-l,n2,-l;ai, 02,03) 



da2 
+Li2fl - 



1 






. 6 


ai + asX 


«! + 02/ 



2 



«! + a2 



+ Li2 1- 



02 + as 
ai + as 



(B5) 

^2 + ttS/^-l 

Other examples together with recursion relations to calculate derivatives may be found in 
Refs. |l6l. Il7|. Some three-electron integrals G can be expressed in terms of T. It is when 
any of argument is equal to zero. Complete list of all cases is: 



G(0, 1,1, 1,1,1) 
^(1, 0,1, 1,1,1) 
^(1, 1,0,1,1,1) 
^(1, 1,1, 0,1,1) 
G(l, 1,1, 1,0,1) 
G(l, 1,1, 1,1,0) 



r(-l, 0, -1; W2 + W3, Wi, U2 + Ms), 
r(-l, 0, -1; Wi + W3, W2, Ui + M3), 
r(-l, 0, -1; Wi + W2, Ws, Ui + M2), 
r(-l, 0, -1; tU2 + Ms, Ui,W3 + U2), 
r(-l, 0, -1; wi + Ms, M2; W3 + Ml), 
r(-l, 0, -1; wi + M2, Ms, m;2 + Mi). 



(B6) 



APPENDIX C: SPECIAL FUNCTIONS 



The complete elliptic integral of the first and second kind, K and E respectively, are 
defined according to [3| as 



K{m) = t dtii-er^'^i-me)-^'^ 

Jo 

(m) = ['dtil-f)-'/^l-mtY^. 
Jo 



E 



They are related to a hypergeometric function 



K{m) 
E{m) 

and fulfill the Legendre's relation 



2Fi(l/2,l/2;l;m), 
2Fi(-l/2,l/2;l;m), 



TT 



E{m) K{1 - m) + E{1 - m) K{m) - K{m) K{1 - m) = -. 



Their first derivatives are 



K'{m) 
E'im) 



E{m) K{m) 
2m (1 — m) 2m 
E{m) — K{m) 
2 m 



(CI) 
(C2) 

(C3) 
(C4) 

(C5) 

(C6) 
(C7) 
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Elliptic functions for |m| < 1 can be conveniently calculated numerically as described in 
1^ . For m > 1 one uses the identity 







E{m ±ie) 
and for m < —1 



mE(-) - 



± i K{1 — m 
1 — m 



K 



m 



\mj 



(C8) 

1 -m) -m)], (C9) 



K 



K{m) = ^== 
E(m) = Vl -mE 



m 



m — 1 
m 

m — 1 



(CIO) 

(Cll) 



The Laurent expansion near the singularity m = 1 is 



K{m) = y: 



n=0 



■ 1 




2 - 


.n\ 


(2)„. 





tpin + 1) - tpin + 1/2) - - ln(l - m) 



(1-m)", 



(C12) 



where ip is a. logarithmic derivative of the Euler Gamma function. 
The dilogarithmic function Li2 is defined by 

ln(l - z) 



Taylor expansion around origin 



Li2(^) = E 



(C13) 



(C14) 



i=l 



is convergent for \z\ < 1. Two useful relations 

Li2(-x) + Li2 ' 



Li2(x) + Li2(l - x) 



TT 

"6 



TT^ In^ X 

'~6 2~' 
— In a; ln(l — x) 



(C15) 
(C16) 



are used for simplification of result of integrations in Eqs. pOl IT^ . Further formulas may be 
found in Q] and an efficient numerical evaluation is described among others in Q| • 
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